Non-linear sigma model approach to quantum spin chains 



in 

o 
o 

o 

Q 

in 



Sumathi Rao 

Harish- Chandra Research Institute, Chhatnag Road, 
Jhusi, Allahabad 211019, India 
email address : sumathi@mri.ernet.in 
(Dated: February 2, 2008) 

We introduce and motivate the study of quantum spin chains on a one-dimensional lattice. We 
classify the varieties of methods that have been used to study these models into three categories, - 
a) exact methods to study specific models b) field theories to describe fluctuations about the classical 
ordered phases and c) numerical methods. We then discuss the J1-J2S model in some detail and 
end with a few comments on open problems. 
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I. INTRODUCTION 



We start with the definition of a spin chainQ, Q as a 
spin model on a one-dimensional lattice - e.g., 



H 



Ji 2J Si-Sj 



-J. 



(1) 



Here, i,j represent the sites on a lattice and and the 
notation nn (=< i, j >) stands for nearest neighbour, 
nnn stands for next nearest neighbour and so on. The 
spins are Heisenberg spins satsifying [Sf , S|] = ie ohc S° 
and not classically commuting variables, and hence it is 
a quantum spin chain. We would like to find the ground 
state and excitation spectrum of these models. 

But why are we interested in these models? Spin sys- 
tems as models of magnetic materials have been used for 
many years because there exist large classes of ma- 
terials where the electron stays localised and magnetic 
properties reside in the individual atoms - i.e., one has 
localised moments which can be modelled by the spins. 

But more specifically, there are several reasons for 
studying one-dimensional quantum spin chain. The first 
is simply that there really exist materials that behave like 
one-dimensional antiferromagnets0, 0. CsNiCl^ is one 
of them, because the ratio between the infra-chain cou- 
pling and inter-chain coupling in this material is 0.018. 
Another compound which is even more markedly one- 
dimensional is NENP {Ni(C 2 H 8 N 2 )2(N02)ClO 4 ) where 
the ratio is of the order 10~ 4 . In both these materials, a 
gap in the excitation spectrum was found although trans- 
lational symmetry remained unbroken. This was an ex- 
perimental verification of a conjecture by HaldaneP,|3,Q 
that S = 1 Heisenberg antiferromagnets should have a 
gap in the spectrum (unlike S = 1/2) and would not 
break translational symmetry (unlike dimers). More re- 
cently, even more exotic compounds which are quasi- 
one-dimensional and can be modelled by unusual spin 
chains (sawtooth spin chains) with missing bonds viz, 
H = J£i&i • S i+1 + J/2£i(l + (-l) l )S, • S 4+2 have 
been found 

The second reason is that there exist exact solutions of 
some toy models, which can then be used as a check or 



testing ground for new anlytical or numerical methods. 
Finally, quantum anti-ferromagnets in higher dimensions 
have become particularly prominent in the last few years 
in the context of high T c superconductors. It is hoped 
that some of the methods to solve quantum spin chains 
may have generalisation to higher dimensions. 



II. VARIETIES OF APPROACHES TO SOLVE 
QUANTUM SPIN CHAINS 

In this section, we will discuss the various methods 
that have been used to 'solve' models of quantum spin 
chains. 



1. Spin- wave theory 

In higher dimensions, the standard way to proceed is 
to start with the classical ground state and then use 
spin-wave theory. We first try to apply that method to 
the onc-dimcnsional spin models here. Let us start with 
the simplest spin-chain, the Heisenberg antiferromagnet 
(HAFM), 



(2) 



Here, i runs over the sites on the one-dimensional lattice. 
If the spins were classical vectors, then 



H 



JS 2 V cos 



(3) 



which is obviously minimum when cos(#i — = 

-1 fa - e i+1 ) = 7T. 

Hence, the classical ground state (Neel state) is given 

by 



>= I U, > 



(4) 



Note that this is not an eigenstate of the Hamiltonian, 
because terms in the Hamiltonian flips nearest neighbour 
spins. However, for very large spins 

[SI Sj] = e a6c S^ = O(S) « 0(S 2 ). (5) 
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Hence, in the limit S — > oo, the Neel state must be 
the ground state. By perturbing about the Neel state, 
we can get the results for large but finite spin. This 
perturbation theory is called the spin-wave theory and 
is done using the Holstein-Primakoff transformation^, 
which is given by 

S? = S-4 ai , S? = -S + b\bi 

S+ = V25(l-||i)V^ , Sf = V2Sbl(l - M*)V» 

ST = V2Sat(l - , S r = V2S(l- b ^/\ 

(6) 

for the A and i? sub-lattices, which are denoted as i E A 
when i is even and i £ B when i is odd or vice-versa. 
We can easily check that the spins satisfy the spin alge- 
bra when the a,i,bi and their conjugates satisfy bosonic 
commutation relations. Note that in the A sublattice, 
the absence of any bosons in a state implies that it has 
the maximum spin and for the B sub-lattice, the ab- 
sence of any bosonic excitation implies minimum spin. 
In the large S limit, the awkward square-root term can 
be dropped and the spin raising and lowering operators 
can be approximated merely as 

Sf -> V2Sa\ , S~ -> V2a* 

S+ -» V2Sb\ , S- -> V2b t (7) 

on the A and B sub-lattices. In fact, we can develop 
a systematic 1/S expansion by expanding the square- 
root term, with the above terms being the first in the 
expansion. But in this review, we will stop with the first 
term. Next, we write the Hamiltonian in terms of these 
bosons (using the above approximation) as 

H = J J2 IS 2 + S ( a hi + b% + a t b, + alb])} (8) 

<»,i> 

After going to momentum space and performing a 
Boguliobov transformation, we get 

H= E E k {c{c k + d\d k ) (9) 

keRBZ 

with^fc = 2JSsin|fc|. Ask — ► 0, E k — ► 2JS\k\, which 
implies that the c and d bosons, which are the spin- wave 
modes, are massless and relativistic modes with spin- 
wave velocity given by v s — 2JS. This, in fact, gives 
us a clue that a relativistic field theory description of the 
spin-wave modes might be possible. 

We can also understand more physically why there are 
two massless spin-wave modes. The Neel state breaks 
the 5*0(3) symmetry of the spin variables down to SO(2) 
(rotations about the S z axis). The spin- waves are the 
Goldstone modes of this spontaneous symmetry break 
down. ( Choosing a direction for the Neel state (ground 



state) spontaneously breaks the 5*0(3) spin symmetry ot 
the Hamiltonian down to SO (2)). 

Spin-wave theory works quite well for three dimen- 
sional magnets, but in low dimensions, spin-wave the- 
ory has problems due to quantum fluctuations. Let us 
calculate the reduction in the sub-lattice magnetisation 
due to quantum fluctuations ( in arbitrary dimensions). 
This can be done by computing the expectation value of 
<5f>. 

< Sf >=< S - a\ ai >=S- a W > ( 10 ) 

k 

In terms of the spin-wave modes, this can be rewritten 
as 

< >= 5 - ^[|wki 2 < 4.c k > +<k < 44 > 

k 

+ K| 2 < > +Uk|«k < c krfk > 

+ Kl 2 ] • (ii) 

All the expectation values are zero in the ground state 
and we are left with 

s,* = s-£ K | 2 ~s-/<g4 (12) 

which is linearly divergent in one dimension and logarith- 
mically divergent in two dimensions. 

Hence, in one dimension, the Neel state is always 
destabilised by quantum corrections. This is just a man- 
ifestation of the familiar result that there is no long 
range order in one dimension (Mermin- Wagner theorem) 
or equivalently, that there is no spontaneous symme- 
try breakdown in 1+1 dimensions (Coleman's theorem). 
Both these theorems are a consequence of the infra-red 
divergences in the theory 

Other methods used in higher dimensions arc fcrmionic 
and bosonic mean field theories. By substituting = 
V'loVH or S, = atoi@ or = ie abc a\ b a ic ^ in the 
Hamiltonian, we get four fermion or four boson terms 
which can then be treated through appropriate mean 
field ansatze. But in one dimension, fluctuations beyond 
the mean field theory turn out to be infra-red divergent. 
Hence, specifically in one dimension, other methods are 
needed. We can divide them roughly into three cate- 
gories. The first one involves the exact solution of some 
model Hamiltonians by some ansatz wave-functions. For 
example 

• Heisenberg AFM for 5=1/2 

The Heisenberg AFM for 5 = 1 /2 in one dimension 
has been solved using Bethe ansatz [ill H^| . The 
solution is hard to write down, but it is known that 
the ground state is unique and that there is no gap. 
Correlation functions fall off algebraically. 
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• S=l model 

The Hamiltonian is given by 

H = • S i+1 - ]T(S, • S. i+1 ) 2 (13) 

i i 

For S = 1, this has a Bethe ansatz solution, which 
shows that the model has a unique ground state 
with no energy gap. 

• Models with valence bond ground states 

— The Majumdar-Ghosh Hamiltonian is given 

by0 

i i 

For S = 1/2, the ground state is given in terms 
of valence bonds. There are two degenerate 
ground states given by 



— = (tl-H) 

There exists a gap in the spectrum and cor- 
relation functions have an exponential fall-off. 
Translational symmetry is broken. 

— The Hamiltonian for one of the Affleck- 
Kennedy-Lieb-Tasaki (AKLT) models [l^ for 
S = 1 is given by 

H = j]rS. t S m + J/3^(SA +1 ) 2 (15) 

i i 

This has a unique valence bond ground state 
found by considering each S — 1 to be built of 
a symmetrised product of 2 S — 1/2's. 




/ = symmetrisation 

The ground state is formed by symmetrizing 
after forming the singlets. Here, again, it was 
found that there exists a gap in the spectrum. 

Besides all these explicit exact solutions of specific 
models, there is another exact statement that can be 
proven in general. That is the Lieb-Schultz-Mattis LSM 
theorem [l6|. This theorem proves that the 1/2 integer 
spin chain either has massless excitations or degener- 
ate ground states corresponding to spontaneously broken 
parity. 



To prove this, let us start with a chain of length L 
obeying periodic boundary conditions. Let us call its 
ground state |i/>o > and assume that this state is rota- 
tionally invariant and an even eigenstate of parity. Now 
construct a new state >— U\ipo > where 

Cr = e (W£)E}=-,Ci+Qs; ) (16) 

i.e., every site from —I to +1 is rotated about the z axis 
through angles in /I, 2in /I, . . .2iir, where I is some num- 
ber of 0(L). First, we have to show that > is degen- 
erate with > hi the L — > oo limit. To do that, we 
compute 

<^\H- Eo\ih >=< MU\H - EoWo > (17) 

where H\ipQ >= Egipo >■ Now using the commutation 
relations of the spins, we can show that 

2 Ttt 2 

< falH - Eolfa >= -^e (2Z + 2) (18) 

where J is the coupling constant of the spins and eo = 
Eq/L. The point to note here is that the R.H.S. is of 0(1) 
and goes to zero as I — * oo. This shows that for an infinite 
chain, \ipo > and \tpi > are degenerate. There is still a 
possibility that asymptotically \ipo >— ► >, so that we 
have only one state. But to disprove that, let us look at 
the behaviour of \ipi > under parity. Under parity, 5? — > 
iSlj and under rotation about the y-axis through tt, Sf — > 
— Sf. Note that both parity and rotation about the y-axis 
through tt, are symmetries of the Hamiltonian. Hence, 
under a combined action of both these symmetries, Sf — > 
-SI,. So 

U = e (iir/£)i:J=-,C7+!)flJ ->u = e (- 27ri )£j = _ ( s? (19) 

Hence, the state >= U\ip >, under a com- 

bined symmetry operation of parity and rotation, 
goes to Ue ( - 2ni ^=-' s i\i/j >= e ^ 2m) ^=-> s i\ip! >. 
But since Ylj=-i &j = + l)^' we see that 

e (w/i) Ej = _,(j+i)SJ = _i if the spin s is odd and is equal 
to +1 if the spin S is even. Hence, for 1/2 integer spins, 
the state > has odd parity and is distinct from \4>o >• 
In fact, < ip \il>i >= 0. Hence, for 1/2 odd integer spins, 
we have proven that as L — > oo, there exists a state \tpi > 
distinct from \tpo >, but degenerate with \ipg >. Hence, 
either there exists a massless excitation with odd parity, 
or if there is a gap, then there is a degeneracy in the 
spectrum. This result is the LSM theorem. The Bethe 
ansatz solution for the Heisenberg AFM with massless ex- 
citations falls in the first class and the Majumdar-Ghosh 
model with two degenerate ground states and massive 
excitations falls in the second class. 

2. Field theory treatment of fluctuations 

The idea here is to derive a low energy continuum 
limit of spin models, keeping only the lowest derivative 
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terms |1|, |6|. We shall first derive the field theory in de- 
tail for the Heisenberg AFM, and then briefly discuss 
how it is done for other general models, including the 
Majumdar-Ghosh model|l3j. 

For the Heisenberg AFM, we start by defining two 
fields 

&2i — &2i+l 



where 



8tt 



(26) 



hi+1/2 = <> x 2i+l/2 = 



25 

&2i + &2i+l 

2a 



(20) 



Here, a is the lattice spacing and the fields are defined at 
a point X2i+i/2 between the sites 2i and 2i + 1 where the 
spins are defined. So the pair of spin variables are now 
replaced by the pair of fields <p an d I- The commutation 
relations for the spins imply that (f>(x) and l(x) behave 
like a scalar field and angular momentum field respec- 
tively. We can also check that $ 2 = 1 + 1/ S-a 2 l 2 / S 2 ~ 1 
in the large S limit. Hence, <j> is a constrained field. 

To derive an effective field theory, we write the Hamil- 
tonian as 

H = J£)E S21 ' S2i+1 + S2 ^ 1 ' S2l] ' (21) 

2i 2i 

then write the spins in terms of the fields and then Taylor 
expand the fields. After doing a lot of algebra, we find 
that 



is an integer ( in Euclidean space). Hence, in the parti- 
tion function, Z = J T>(f>e~ s , e l6 ® = e 2mS< ^ is periodic in 
5. 5 = is equivalent to all 5 = integers and 5=1/2 
is equivalent to all 5=1/2 integers. Also, we note that 
for integer spins, the topological term can be dropped 
because 27TlS< 3 =1 for all configurations, but for half in- 
teger spins, it is either +1 or —1 depending on value of 
Q. Thus the Hopf term plays an important role for half- 
integer spins. This was what in fact, led to the famous 
Haldane conjecture that the HAFM for integer spins has 
a gapped spectrum and is massless for half-integer spins. 

From the field theory mapping, in fact, it is easy to see 
that integer spins models have a gap, but it is more non- 
trivial to show that half- integer spin models are gapless. 
Let us start with a semi-classical analysis of the integer 
spin models. Semiclassically, we assume that the 5(3(3) 
symmetry of the Lagrangian is spontaneously broken to 
U{1) ~ 50(2) by the Neel state or vacuum state given by 
4> — (0, 0, 1). Fluctuations about this state are described 
by 



2 ,(l-0 2 -0 2 )72) 



L,02,l) (27) 



f - 5 S 2 <j)' 2 

H = 2Ja / dx[(l + -^4>') 2 H ^ — ] (22) to linear order in fluctuations. Hence, the Lagrangian 



dx 



and X^2i(2 a ) — J dx. We now introduce L 



where q 

the spin-wave volocity v s = 2JaS and also the coupling 
constants g 2 = 2/5 and 9 = 2irS. This allows is to 
rewrite the Hamiltonian density as 



1 

2^ 



2 5^0! 



1 

2fJ 



2^02^02 (28) 



2 An g z 



(23) 



which, with some more algebra can be shown to be de- 
rived from the Lagrangian density given by 



d^cj> + — e» v <p ■ x d v <i> (24) 



with 4> 2 = 1. Note that we have already taken the large 
5 limit. This is necessary not only to have 2 = 1, but 
also to justify the Taylor expansion. By keeping terms 
only upto second order in derivatives, we are assuming 
that the deviations from the equilibrium positions of the 
spins are small, which is justified only in the large 5 
limit. With these assumptions, we find that the spin- 
wave modes or fluctuations in the HAFM are described 
by an 0(3) non-linear sigma model (NLSM) with a Hopf 
term (the term proportional to 9). 

The Hopf term is a total derivative, but its integral is 
an integer. Hence, the action 



is just the Lagrangian of two free bosons. This is the same 
as the result that was obtained using spin-wave theory. 

But using the field theory, we can do a lot better. 
Firstly, we can use renormalisation group (RG) to go be- 
yond naive perturbation theory, i.e., we can compute the 
/3-function. Since the manifold here ( of values taken by 
the fields {4>i,4>2,4>?,)) is a sphere, we can use geometric 
methods to compute the RG equation and we find that 



/V) 

9 2 eS (L) 



dg 2 



9 



d\nL/a 2n 



(1 - G#nL/o)/27r) 



(29) 



5 = 



dtdxL — — - I d 2, j ' 



+ it 



(25) 



where g 2 is the microscopic coupling that was derived at 
length scale L — a to be 2/5. From this, it is clear that 
the coupling constant blows up when (g 2 lnL/a)/27r = 1 
which implies L/a — e 27r ^ 9 — e* s . Thus, as a function 
of g 2 , we expect a phase transition to the strong coupling 
regime, where the earlier perturbative result of two mass- 
less bosons is no longer valid. Since the length scale is of 
0(e ), masses of order 0(e~ vS ) are expected i.e., one 
expects to flow to a strong coupling regime, where there 
is a gap of 0(e~ ) to excitations. 
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One can also substantiate this by solving the field the- 
ory in the large N limit, i.e., by extending the 0(3) 
NLSM to O (TV) [HQ, with a Lagrangian 



(30) 



4> 2 N = 1. In other words, 



with 4> 2 = 4>\ + 4>l + <f> 2 + . 
instead of having just the usual spin variables with three 
components, we have extended it to N components. This 
can also be thought of as taking the number of dimensions 
in which the spin moves to be N. In the limit of large 
N, it is actually possible to compute the path integral 
explicitly and obtain the mass generated and we find that 



(31) 



for each of the N bosons, where A is an ultra-violet cutoff. 
As N — ► oo, S — * oo, but A — * oo as well, so as to keep m 
fixed. Higher order corrections will go as 0(1/N). Hav- 
ing obtained this result for large N , we now bravely set 
N = 3 ( assuming corrections will be small) and conclude 
that the integer spin HAFM has an excitation spectrum 
consisting of a triplet of massive bosons with masses of 
the order of e~ vS . 

All of this was for integer spins. Now what about 1/2 
integer spins? Here, the field theory includes the non- 
trivial Hopf term and is quite difficult to solve. How- 
ever, Affleck 2] has mapped the model to a k = 1 Wess- 
Zumino-Witten (WZW) model and by studying its sym- 
metries, he has argued that the 6 = it case is massless. 
This difference between the integer and half-integer spins 
was the big contribution of field theories in spin models. 

Similar mappings have also been used to write down 
field theories of other models, such as the Majumdar- 
Ghosh model and its generalisations^! For instance, 
for the MG model for arbitary spi ns, we can write down 
an S0(3)l X S0(2)r field theory [l5j by introducing an 
50(3) group valued R field as follows - 




R = 



In terms of the R field, the Lagrangian is given by 
1 



L 



T tr{R R) 



tr(R T R'I 2 ) 



4cg 2 v ; 2g 2 



(32) 



(33) 



with g 2 = V6/S and c = JSa a/27/8 and I 2 being a 
diagonal 3x3 matrix with diagonal entries (1, 1,0) and 
all other entries zero. Here, R denotes the time derviative 
of the matrix- valed field R and R', its space derivative. 
The fields <f>i are related to the spins as 

S3i_l — S3i + i 



(01 ) 3* 

(#2)3* 



V3S 

831-1 -I- S 3i+1 — 2S 3i 



V35 

3)3i = (<£l)3i X (^2)3i 



(34) 



Note that the field theory has no topological term. This 
is not unexpected, because here the manifold of the fields 
is 50(3) and n 2 (50(3)) = 0, whereas for the HAFM, 
the manifold was S 2 and H 2 (S 2 ) = Z. So at least naively, 
no difference is expected for integer and half-integer spin 
models. Also, note that the global symmetry of the ac- 
tion is SO(3)l x SO(2)r, which means that the effective 
action at any length scale can be written as 



.L )tr(R ^ ) + ( _l__I s 



2 )trR T R'I 2 



with the microscopically derived Lagrangian having g\ — 
92 = 93 = 2^1 = g 1 = \/Q/S. But these values change as 
we go to larger length scales in accordance with the RG 
equations or /3-functions given by 

2 / 1 1 M 

■7 : r + 5i 93{— 2 ) 

(3i 94 + 3233) 9i 32 





gi 


9\ ^ 


2ir 








9i 


9l = 


2 k 








gi 


3 3 2 = 


2ir 








9t 


9l - 


2ir 



9l 



9i 93 (-5 
3i 

r,1 n r. 



-) 2 + 23i 3s( 

32 



<i 



1 



— 7 1 s+9l93{ 2 

9i (3i34 + 3233) 33 



9l 



9s 9i (-2 

33 



Si' 



231 33 ( — 
34 



9l 



(35) 



We integrated these equations numerically^] and found 
that the length scale where strong coupling takes over is 
Q = L/a = e 5 - 76S , which is of the same order as e* s that 
we had found for the HAFM. Moreover, we found that 
the flow is such that 31/32 and 33/34 flow to unity, so 
that the symmetry gets enhanced to 50(3)^ x 50(3)^, 
and Lorentz invariance is restored. Thus, the Majumdar- 
Ghosh model for arbitrary values of the spin flows to a 
disordered phase. We shall come back to this analysis in 
the last section where we study a general dimerised and 
frustrated model. 

3. Numerical methods 

The third method that has been used to study spin 
chains is through numerical computation. Here, I shall 
only quote various results. 

• Exact diagonalisation of small systems 
The frustrated Heisenberg antiferromagnet mod- 
eled by 



(36) 



has been studied for S = 1/2 to upto 20 sites and 
it was found that the critical value of a for which 
a gap opens upin the spectrum is give by a c = 
.2411 ± .00010- This is the point at which the 
fluid-dimcr transition takes place. 
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• Density-matrix renormalisation group 

(DMRG) [fj| 

This is recent method which has gained ground 
and is remarkably accurate. The idea is to combine 
exact diagonalisation methods with the idea of 
renormalisation group. So a small system is first 
diagonalised exactly and then the system size is 
increased by adding two spins at a time on cither 
side. This is done repeatedly using RG ideas. 
For the same model as above, DMRG also finds 
a c — .241. DMRG was also used to study a more 
general model involving bond alternation 20J. 




0.25 j 1.0 

FIG. 1: Semi-classical phase diagram of the J 1 — J2 — S model 



III. FRUSTRATED AND DIMERISED AFM SPIN 
CHAIN 



The idea is to study the Ji-J 2 -S model given by 

H = Jx ]T[1 + (-iys}s l s l+ i s A+2] ( 37 ) 

i i 

in detail^?! Classically, the ground state is a coplanar 
configuration of the spins with energy per spin 



Jl 



E o = 5^[y(l+<J) cos9 1 + A-{l-6)cos9 2 +J 2 cos(9 1 +9 2 )} 

(38) 

Minimising this energy with respect to 9i gives three 
phases 

• Neel phase 



nnnri--. (39) 

This is stable for (1 - 6 2 ) > 4J 2 /Ji. 

• Spiral phase 

Here, the angles between neighbouring spins alter- 



\ f 



nate between 9\ and 9 2 where 

1 A- 5 2 8 4J 2 

cos 9\ 



and cos 9\ 



l + S l 4J 2 /Ji 1 + S 2 Ji 
1 A-S 2 S 4J 2l 



(40) 



1-(T4J 2 /Ji l-<5 2 Ji 

This phase is stable for IS 2 < 4J 2 /Ji < (l-8 2 )/5. 

• Collinear phase 
This phase can be thought of as a special case of 
the spiral phase where 9± = ir and 9 2 = 0. It can 
be denoted as 



mmu 



(41) 



This phase needs both frustration and dimerisation 
and is stable for (1 - 5 2 )/S < AJ 2 /Ji. 



These three phases in the classical phase diagram are 
depicted in Fig.(l). 

We can study fluctuations about the classical ground 
state as described earlier. In the Neel phase, there are 
two modes with equal velocity and the Fourier transform 
of the spin-spin correlation function S(q) is peaked at q = 
7r. In the spiral phase, we have three modes, two of them 
with equal velocity describe out-of-plane fluctuations and 
the third one with a higher velocity describes in-plane 
fluctuations. S(q) is peaked at ir/2 < q < it. In the 
collinear phase, once again, there are two modes with 
equal velocity, but here S(q) is peaked at q — Tt/2. But 
as we have already seen earlier, we do not expect spin- 
wave theory to be accurate in one dimension because, 
there is no long-range order, no spontaneous symmetry 
breakdown and no Goldstone modes in one dimension. 

Next, what do we know about the model exactly? For 
J 2 = S = 0, the model is just the HAFM and the solu- 
tion for S — 1/2 is a unique ground state with no excita- 
tions. For J 2 = Ji/2 and 5 — 0, which is the MG model, 
the solution for S = 1/2 is the doubly degenerate valence 
bond state, with massive excitations. In fact, this state 
turns out to be the ground state even with dimerisation 
along the line 2J 2 + 5 = J\. 

Now, let us study the field theory model for the fluc- 
tuations in the three classical phases. 

• In the Neel phase, even with J 2 and 5, the map- 
ping is to an 0(3) NLSM, with the Hopf term as 
given in Eq.(j2U- The only difference is that now 
c = 2JiaSyi -S 2 -4J 2 /Ji, g 2 = 2/(5(1 - S 2 ~ 
4J 2 /Ji)) and 9 = 2wS(l-S). We expect the theory 
to have a mass gap in general and to be massless 
only when 9 = 2^5(1 — 6) = tt. Note that a topo- 
logical term is present to distinguish different spins, 
but spin is not really a continuous variable. So for 
each spin, integer or half-integer, there are specific 
values of 5 which can be chosen to get massless 
points. 

• For the spiral phase also, the field theory still turn 
out to be the SO(3)l x SO(2)r invariant, but with 
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a Lagrangian given by 
1 



1.0 



L = 



tr(R T RP ) - Tr-yHR T R'Px) (42) 



4cg 2 v u/ 2g 2 



where Pq and P\ are diagonal matri- 
ces with the diagonal elements given by 
P = (l/2 3 2 ,l/2 5 2 ,l/ 5 2 i/2«,2) and 
Pi = (l/2 5 2 ,l/2. 9 2 ,l/ 3 2 - l/2 ff 2 ) respectively. 
The RG equations are the same as the ones given 
in Ea. (|35|l . However, the initial microscopic values 
of the coupling constants are different now and are 
given by 



fll 



9l 



4 J 2 + Ji 
4 J 2 - Ji 



and 



(1 - Ji/2J 2 ) 2 ] . (43) 

As before, the RG equations can be integrated nu- 
merically with these initial conditions and it can 
be shown that the theory flows once again to an 
SO(3)l x SO (3) a Lorentz invariant field theory. 
The interesting point is that this theory turns out 
to be an exactly solved model j^. The low en- 
ergy spectrum consists of a massless spin 1/2 dou- 
blet. Hence, in the spiral phase ( which requires 
sufficiently large frustration and dimerisation) , we 
can make the prediction that both integer and half- 
integer spin models should have massive spin 1/2 
excitations. The long wavelength excitations are 
expected to be 'two-particle' excitations, the spin 
triplet and the spin singlet excitations. 

Although there is no topological term in the La- 
grangian, we claim that there does exist a differ- 
ence between integer and half-integer spins in this 
phase. Tunneling between soliton sectors can lead 
to a unique ground state for integer spins, but this 
is not possible for 1/2 integer spins, which have a 
doubly degenerate ground state, in accordance with 
the LSM theorem. 

• Finally, we can write down the field theory for the 
collinear phase as well. Here again, the field theory 
turns out to be an 0(3) NLSM, but without the 
Hopf term. This means that the phase is always 
gapped both for integer and non- integer spins. 

We generally expect these field theories to be valid in 
the large S limit, but for small values of S such as S — 
1/2 and S = 1, the above analysis is only indicative and 
numerical studies are needed to get thephase diagram 
accurately. These have been obtained [la, |2fJ and we only 
reproduce the phase diagrams here - 

As can be seen by comparing these diagrams, with the 
classical phase diagram in Fig. (1), the qualitative pic- 
ture is reproduced for the spin 1/2, but for spin 5 = 1, 
there are many new unexpected features in the S = 1 
phase diagram obtained numerically. 
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FIG. 3: S — 1 phase diagram of the Ji — J2 — S model 



IV. CONCLUSION 

In this paper, we have given an overview of the field 
of quantum spin chains, with emphasis on the non-linear 
sigma model mapping. To recapitulate, quantum spin 
chains are spin models on a one-dimensional lattice. For 
parity invariant systems, the Lieb-Schulz-Mattis theo- 
rem says that for half-integer spin models, the ground 
state is either doubly degenerate, or the spectrum con- 
tains a massless mode. Using the NLSM mapping, we 
demonstrated that the difference between 1/2 integer 
spin chains and integer spin chains was caused by the 
existence of a topological Hopf term in the Lagrangian. 
The presence of this term for 1/2 integer chains led to 
a gapless spectrum, whereas integer spin chains which 
did not have the Hopf term were gapped. For more 
general models, such as spin chains with dimerisation 
and/or frustration, the NLSM approach can only give a 
qualitative understanding. For instance, the mapping of 
the Majumdar- Ghosh model (more generally, the spiral 
phase of a frustrated and dimerised spin chain to the RG 
fixed point Lagrangian of an 50(3)^ x SO(3)l model 
leads to the prediction that the low energy spectrum 
consists of a massive spin 1/2 doublet. But for low val- 
ues of S, such as 1/2 and 1, often numerical methods 
are needed to get better results, as seen in the explicit 
phase diagrams for the spin 1/2 and spin 1 frustrated and 
dimerised models . 
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One of the important issues in this field is to get a 
proper understanding of the Haldane gap. Usually, a 
gap is formed when some symmetry is broken. So we 
need a symmetry that exists for half-integer spins and is 
broken by all integer spins. Since the distinction between 
the integer and half-integer spins occurs because of the 
topological Hopf term, it is expected that the order pa- 
rameter characterising the massive and massless phases 
is also topological in nature. A claim is that there exists 
a hidden Z2 X Zi symmetry in the S = 1 model, which 
when broken leads to the gapped Haldane phase. But 
this phenomenon is not well-understood. 
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